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ABSTRACT 


The excitation function for the Shee x) pro reaction has been 
experimentally determined in the proton energy range from 1420 kev 
through 2160 kev. Unreported resonances in this reaction have been 
identified at 1975, 2033, 2075, and 2117 kev with estimated errors 
of + 4 kev. Gamma spectra at these resonances are obscured by 
contaminant radiation but suggest complex decay schemes. Double 
angular correlation measurements have been made at the 14/70-kev re- 
sonance indicating a resonance level spin of 2+, or less probably, 
2-. The level at 4.50 Mev is shown to probably have a spin of 2+, 


or less probably, 3+ rather than 0+ as has been predicted. 
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I. INTRODUCTION 


As a result of successive refinements and extensions of the collec- 
tive model of nuclear structure, it has been possible to extend its ap- 
plicability into the region of intermediate mass nuclei (20€A40). By 
considering the outer nucleons to be strongly coupled to the nuclear core, 


excited states in Hila eo me and ee have been interpreted as col- 


lective rotational excited states with considerable success. More re- 
cently states in Pa and T° have been identified with collective vi- 
brational states of a weakly coupled collective nodes The experimental 
collection of information on the states of intermediate mass nuclei serves 
to verify the theoretical models, to determine the parameters which enter 
into descriptions of the models, and to permit the extensions and refine- 
ments to theory necessary for a more adequate description of nuclei. Of 
particular relevance to the present investigation is the prediction based 
on the weakly coupled collective model applied to A = 30 nuclei that the 
4.50-Mev level of ae has a spin-parity of ee In the present work the 
experimental determination of this value has been the object of double 
angular correlation measurements performed on the gamma rays observed at 
the 1470-kev resonance in the ST OE reaction. 

Determinations of resonance energies in the Staa(p, Chee reaction 
have been carried out by numerous investigators since Cre The pro- 
ton energies which have been employed in these investigations extend over 
the range 300 to 1950 kev. 

In addition to detecting resonances and measuring the energies at 
which they occur, the allied problem of analyzing the associated gamma-ray 
spectra has been the object of an increasingly large amount of research 
commensurate with improvements in gamma-ray spectroscopy techniques and 
instrumentation. 

Investigations of the gamma rays which arise from deexcitation of 
the excited states of a were begun in 1954 by Endt, Kluyver, and Van 
der rea” They studied the Si (peux reaction at resonances oc- 
curring at proton energies of 414 kev and 326 kev and assigned spin, 
parity, and isotopic spin to the ground state and the first excited state 
of p>? premised on gamma-transition probabilities, shell model consider- 


ations and the isotopic spin multiplet formalism. In 1956 Broude et alee 
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confirmed the assignments to the first excited level by measuring the 
angular distribution of the gamma-rays at the 414-kev resonance. They 
also detected five new states in Po by examining gamma radiation 
arising from the same reaction at resonances occurring at proton enere 
gies of 696 kev and 737 kev. 

Investigations by Lee and eee” and by Endt and pees who 
studied the s?*(d,a)P°° reaction, showed the first excited state to be 
in reality a close doublet. The latter investigators were able to accu- 
rately determine the energies of 30 excited levels of peo up to excitation 
energies of 5.8 Mev by magnetic analysis of the alpha-particle groups from 
Chis teaction. 

In 1958 Van der Leun and pude= reinvestigated the four aforemen- 
tioned resonances, studying the decay schemes by analysis of gamma-ray 
spectra, coincidence measurements, and angular distribution measurements. 
From this work came decay schemes, branching ratios, and assignment of 
spin-parities and isotopic spins to the four resonance levels and to five 
other excited levels of pr? including the two states of the above men- 
tioned doublet. During 1962 the gamma-ray spectra associated with the 
proton resonances at 1330 peas and the resonances at 13/5 kev and 1505 
kev were ieee eae Proposed decay schemes were advanced for these 
resonances and spin-parity assignments made to three levels of ee One 
of the more comprehensive investigations of the sia ye reaction 
is that of Baart et aio who studied six resonances and proposed decay 
schemes for them. Firm spin-parities were assigned to six low-lying 
levels and tentative assignments were made to four more levels. The 
most recent investigation of the reaction is that by Meera As a re= 
sult of this work decay schemes were proposed for resonances at 1302, 
1470, 1505, 1640, 1686, 1748, and 1772 kev, and spin-parities were de- 
termined for three resonance levels and two intermediate levels. The 
present knowledge of the quantum states of the Lower Bae levels re- 
sulting from the above investigations is summarized in Fig. 1. Two come- 
pilations which summarize important experimental results of the above 
mentioned investigations are the Nuclear Data shesuees and the compi- 
lation of information on the energy levels of nuclei with Z = 11 to 


Z = 20 by Endt and Van der raging 
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The purpose of the present work has been to investigate resonances 
in the ee a ye reaction in the proton energy range from 1950 kev 
to 2150 kev which has not previously been reported on, and to investigate 
the gamma-ray spectrum associated with the known resonance occurring at 
1470 kev. In particular, the angular distributions of gamma rays at the 
latter resonance have been measured in a double angular correlation ex- 
periment to obtain information on the 4.50-Mev level in oe Additionally, 
some information has been obtained from examination of spectra taken at the 
higher resonances although contamination by ie hampered this work. The 
following sections of this paper describe the experimental equipment, pro- 


cedure, and data analysis employed, and the results which have been ob- 


tained from the investigation. 





II. DESCRIPTION OF EQUIPMENT 


1. Experimental Layout 

Fig. 2 indicates the general arrangement of equipment used in this 
experiment. Fig. 3 is a photograph of the target area showing the cor- 
relation table tilted into a vertical plane. 

2. Accelerator 

The 2-Mev Van de Graaff horizontal electrostatic accelerator of the 
U.S. Naval Postgraduate School, manufactered by the High Voltage Engi- 
neering Corporation, was used in this work as a source of protons. The 
radio frequency ion source is of the type described by o.°° A proton 
spin-resonance fluxmeter for accurately measuring the analyzing magnet's 
field was installed furing the course of the work. An error signal from 
the fluxmeter was used in a feed back loop to the power supply of the 
beam analyzing magnet to hold the field strength automatically at a set 
value. This feature simplified the maintenance of a given resonance 
energy for long bombardment times without "slipping off". 

A current integrator of the vacuum tube electrometer type was used 
to measure beam current. High stability plastic insulated capacitors 
were used to collect the charge. The voltage across the capacitor in 
use was applied to the grid of a 954 acorn tube, whose output fired a 
Schmidt trigger circuit at a preselected voltage level. The trigger 
output actuated other circuitry to discharge the capacitor and advance 
a Veeder-Root cycle counter after each cycle. 

3. Target Assemblies 

Two types of targets and target mounts were developed and used in 
this experiment. These are shown in Fig. 4. For resonance and coinci- 
dence measurements, which were made in the plane perpendicular to the 
beam, the target backing consisted of a one-half-insdiameter flanged 
copper hemisphere held against a flat copper boss on the mount by a 
threaded retaining ring. This assembly was designed with the beam di- 
rection as a symmetry axis so that it can be used for triple correlation 
experiments in the perpendicular plane with no perturbations from target 
geometry. For the double correlations, measured in the horizontal plane 
containing the beam, target backings were 0.046 cm thick copper strips 


clamped in a copper mount. Both types of backings and mounts were vacuum 
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plated with a layer of gold thick enough to stop incident protons. For 
both geometries the target mount was fitted inside one end of a one-in. 
diameter tubular chamber, the axis of which coincided with the beam path. 
Both target assemblies were designed so as to keep mass in the target 
vicinity to a minimum. Competing with this requirement was the need to 
cool the target in order to prevent any appreciable target material eva- 
poration. Maximum possible beam intensities were desired so as to take 
advantage of the better counting statistics which accompany high yields. 
Still another requirement of angular correlation measurements is a small 
target spot to approximate a point source of gamma rays. A beam colli- 
mator consisting of two fiveemil tantalum discs separated by 0.4 in. was 
mounted in a brass holder 38 in. from the target. A 0.096-in. diameter 
collimating aperture was followed by a 0.098-in. hole in the second disc 
for the purpose of stopping protons scattered from the edge of the cole 
limating disc. A cylindrical lead shield surrounded the collimator in 
order to attenuate any gamma rays produced here. 

As is described in Chapter III, target material was evaporated di» 
rectly onto the gold-plated copper backings to permit good heat conduction. 
Initially, a hollow copper target backing through which coolant could be 
circulated was tried, but was abandoned for a simpler solid copper back- 
ing cooled by conduction through a finned copper extension leading to an 
external Dewar flask. Liquid air in the Dewar proved to be convenient 
as a coolant at high beam powers and, below about ten watts, target power 
dissipation by forced air cooling of the target extension fins proved 
adequate. 

After unsuccessful trials of Lucite and glass target tubes, aluminum 
was finally adopted for this purpose. An auxiliary target assembly with 
a quartz window target was used to accurately center the proton beam spot. 
The target chamber tube, which was insulated from ground at both ends, 
was maintained at a negative potential of 300 voits relative to ground in 
order to suppress secondary electrons fron the target. in order to pres 
vent recoil protons from exciting aluminum resonances, the target tube 
was lined with one-mil gold foil. 

The hemispherical copper target backings for measurements in the 
plane perpendicular to the beam were needed in some quantity, so were 
"mass produced" in a small hand forge made for this purpose. A surface 


which needed only minor polishing and cleaning before 
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A frequently troublesome feature of (p,% ) investigations is the 
gradual accumulation of carbon on the target. as a result of cracking 
of pump oil vapor by particle bombardment. To reduce this difficulty, 
advantage was taken of an 8-in. glass cross which was mounted between 
the target diffusion pump and the target chamber as shown in Fig. 5. 

A special liquid air cold trap was constructed and installed in this 
cross to serve a dual purpose. First, it functioned as an ordinary 
cold trap with a large surface of copper vanes suspended from the 
coolant reservoir to provide an optically dense path from the dif- 
fusion pump thus permitting the vacuum to be routinely maintained close 
to iC Torr. The second feature of the trap was a lein. diameter 
copper pipe piercing the vanes which the proton beam traversed. At 
high vacuum the proton beam functions effectively as 4 pump, trans-~- 
porting gases and vapor encountered in its path. When the beam passes 
through the restricted cold pipe for 6 in. before entering the target 
chamber, it was reasoned that its drifting retinue of condensable vapor 
contaminants would be largely trapped by the cold inner wall. The re- 
sults from use of the trap support this idea since only minor discolor- 
ration was observed on targets compared with that which occurred when 
the trap was not used, and no increase in bombardment energy was found 
necessary as targets aged. 

4. Shielding 

To reduce the radiation level at the target resulting from brems- 
strahlung and X-rays caused by electrons accelerated toward the high 
voltage accelerator terminal, a considerable amount of lead shielding 
was erected. A 2-in.-thick lead wall, 5 ft. high and 4 ft. wide, was 
placed between the energy control slits and the beam collimator. In 
addition, two lead screens 1/4-in. thick and each 5 ft. by 4 ft. were 
bolted onto a caster supported steel framework constructed for this 
purpose. The mobility of the screens permitted theiz placement to 
minimize the radiation leve! at the target and also allowed easy re- 
moval for access to the target area. Although the above-described 
shielding offered some reduction in natural background radiation, 
(notably the 1.46-Mev Ce gamma ray) it was not feasible to screen 
specifically against background without either increasing scattered 
target radiation at the detectors or utilizing an excessive amount of 


lead. 
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5. Correlation table 

A correlation table shown in Fig. 6 was built for this experiment 
to support the gamma-ray detectors and to permit their accurate positioning 
with respect to the target. A steel base plate mounted on casters supported 
a hollow magnesium column upon which the table was mounted. Screw jacks in 
the base plate were extended when the apparatus was in place to hold the 
table in position and to allow accurate vertical adjustment. The table ite 
self consisted of a movable graduated aluminium disc which was mounted con- 
centrically on a fixed aluminum support plate. Slides mounted on the disc 
and on the support plate held mounts for the two gamma-ray detection as-~= 
semblies. The mounts were positioned in the slides by threaded rods which 
allowed the detectors to be set accurately at distances from 1/2 in. to 
6 in. from the target position. The detection assembly mounted on the disc 
could be positioned in angle with respect to the fixed detector in 12 min. in- 
crements of arc by means of an adjusting gear meshing with gear teeth on 
the periphery of the disc. The entire table assembly was attached to the 
magnesium support column by a hinge so that it could be tilted into either 
a horizontal or vertical plane, permitting measurements to be taken in 
either the horizontal plane containing the proton beam or in the plane per- 
pendicular to the beam. The first named arrangement is designated herein 
as geometry A and the latter as geometry B. 
6. Gamma-ray detectors 

The gammaeray detection assemblies were supported with their center- 
lines 8 3/4 in. above the table. They each consisted of a 3 in. x 3 in. 
cylindrical sodium iodide (thallium activated) crystal integrally mounted 
with a Dumont 6363 photomultiplier tube supplied by the Harshaw Chemical 
Company. In the earlier phases of this work, similar 2 in. x 2 in. crystals 
mounted with Dumont 6292 tubes were used. The resolution of the 3-in. cry- 
stals, measured after the signal had traversed a l2-microsecond delay line 
input to the pulse height analyzer, was 9% for the 0.662-Mev conae gamma ray. 
A 6AK5 cathode follower was constructed into each detector mount to serve 
as a preamplifier. 
7. Electronics 

The electronic equipment and circuitry for processing signals from 
the detectors were planned so as to permit considerable flexibility in exe 
perimental applications. A block diagram of the basic layout is shown in 
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The system is composed of an RCL 128-channel pulse height analyzer 
to record pulse height spectra, a slow-fast-sum coincidence analyzer, and 
scalers for recording various outputs. 

It was found that requirements for the detector output signal by 
several components, including a monitor oscilliscope, tended to overload 
the preamplifiers and to permit unacceptable interaction of subsequent 
Stages with each other through the coupling of their inputs. To avoid this 
problem a bank of six cathode followers was constructed. A bussing arrange- 
ment made it possible to feed the grids of the followers from either one of 
two busses or from an individual input, so that up to six pieces of equip- 
ment could be driven by the same input pulse with no interaction and with 
practically no signal loss. Each follower was provided with an individual 
output gain control. 

The output of either detector as desired was also fed directly into 
the pulse height analyzer through a simple resistance isolation network. 
From the resistance network the signal passed through a terminated 1350-ohm 
delay line to allow time for formation of gating pulses needed for coinci- 
dence spectra. Several methods of injecting a signal into the analyzer 
Were tried inasmuch as this equipment was designed to accept an input only 
from a photomultiplier tube connected directly into the preamplifier as- 
sembly provided. This was not an acceptable practice because there existed 
no convenient means of signal extraction for other purposes, and because 
the preamplifier utilized 4a high-voltage supply which required running the 
photomultiplier cathode at a high negative potential with respect to the 
tube shield with the attendant probability of tube damage. The procedure 
adopted as most successful was to apply the signal from the delay line to 
the anode pin of a l]4=pin tube base which was connected directly into the 
analyzer preamplifier. It was also necessary to incorporate a gate pulse=- 
forming network into the analyzer in order to meet the gating pulse speci- 
fications required by this equipment. 

Gamma-ray spectra were obtained by gating the analyzer with the dise- 
criminator output pulses trom the linear amplifier fed by the same de- 
tector. This procedure excluded a weak precursor pulse which was trans- 
mitted by the delay line inpit, by not gating the analyzer on until a few 
microseconds before arrival of the delayed pulse. No significant loss in 


detector resolution resulted from leaving the delay line in the input for 
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ordinary gamma-ray spectra, and it was found quite convenient to be able | 
to switch from these singles spectra to coincidence spectra without change | 
in analyzer calibration. 

For the observation of coincidence spectra, the analyzer was gated by 
the slow-fast coincidence output. The basic slow-fast coincidence system 
assembled for this experiment has been described emeheree” In this 
System the output from an Atomic Instrument Company slow coincidence unit 
accepting up to four inputs was combined in a diode mixer with the fast 
coincidence unit output to obtain the slow-fast coincidence output signal. 
The fast unit, based on the circuit of Draper and pigts@fer, 2 detects coe 
incidences between amplified pulses from the scintillation detectors with 
a resolution time as low as ten ns. For pulses spanning a large amplitude 
range, however, more dependable operation was obtained when the resolving 
time was increased to about 80 ns. This was due to an observed amplitude 
dependence of the shape of the leading edge of the amplified pulse, to 
which the fast unit is sensitive. 

Several modifications have been made to the coincidence system which 
do not appear in the previous description. The diode mixer unit for come 
bining the slow and fast coincidence outputs has been modified by inserting 
a Schmidt trigger between the adding stage and the output stage in order to 
clip the sum pulse and improve output-pulse characteristics. The fast unit 
has been modified by providing a variable delay line in the input in order 
to permit accurate compensation for differences in time delays between the 
detectors and the coincidence unit. Finally, there has been added to the 
system a provision for obtaining sum-coincidence spectra as first described 
by Booeenboon.”° In this system the pulses from the two detectors are 
carefully matched in amplitude calibration, then summed in a pulse adding 
circuit. These summed pulses are analyzed in a differential discriminator 
set to accept pulses corresponding to the excitation energy of a resonance 
level. This output is then used as a third input to the slow coincidence 
unit so that its output represents coincidence between two completely ab- 
sorbed gamma rays comprising a double cascade from the resonance level. 
This signal is then used to trigger the pulse height analyzer gate. The 
spectrum obtained from pulse height analysis of one detector output thus 
consists only of full energy peaks due to radiations which are members of 
a double cascade. The pulse height adder which has been constructed for 
this purpose employs the circuit used by Kantele and Fink in a somewhat 
different type of Haein The adder output was analyzed by a con- 


ventional linear amplifier and discriminator. 
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The scalers used to record single channel counts had dead times of one 
microsecond and those recording coincidences had dead times of approximately 
eight microseconds. Separate cathode-follower outputs to drive scalers were 
constructed for all coincidence stages so the relatively low scaler input 
impedances would not degrade coincidence pulses. The number of scalers proe= 
vided, which exceeded the requirements for the present experimental measure- 
ments, permitted a continuous check on the performance of the coincidence 
system. 

A double pulser with associated pulse shaping networks and an os= 
cilloscope formed an integral part of the equipment and were used for 
all routine checks and adjustments. An additional slow coincidence unit 
of variable resolving ne which was constructed as part of the system, 
was used to check on the operation of the Atomic Instrument Company unit. 

Control circuit relays were provided so that the pulse height analyzer 
and all scalers could be started and stopped simultaneously by the current 


integrator on-off switch. 
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IIIl. EXPERIMENTAL PROCEDURE 


l. Target preparation 

The target material used in this experiment was isotopically enriched 
silicon in the dioxide form purchased from the Stable Isotopes Division of 
Oak Ridge National Laboratory. Relative percentages of the three stable 
Silicon isotopes were 91.85 ssi, 7.85 igen and 0.30 Si This material 
was evaporated onto prepared target backings in a vacuum evaporation ape- 
paratus under a chamber pressure of 10°? to 10°° EOET 

Beam power of as much as 20 watts was delivered to a target spot area 
of about 6 square millimeters during bombardment. Dissipation of this heat 
by 5-mil tantalum target backings, which had been used in previous work at 
this laboratory, was found insufficient to prevent local overheating of the 
target with an accompanying rapid loss of target material. This fact, and 
the fluorine contamination of the surface of all sheet tantalum examined, 
led to employment of the heavier target backings of copper described in 
Chapter II. While vacuum plating these copper backings with 60,000 K or 
more of gold (99.999 percent purity) to completely stop 2-Mev protons, it 
was found that the plating tended to crack and spall off when the copper 
backings were maintained at ambient temperature during evaporation. This 
effect was caused chiefly by tensile stresses which have been found to be 
as much as 850 Ke in 2000 & gold films deposited on eoeren > This 
trouble was overcome by holding the copper backings at 400°c during evae= 
poration by means of a resistance heater which was constructed for this 
purpose. The heater was made in the form of a 4-in. cylindrical section, 
to which the target backings were clamped. The cylindrical heater section 
was suspended in the vacuum coating chamber so that its axis coincided 
with the evaporant source. Heating the backings in this manner prior to 
evaporating both gold and $i0, also served to remove any residual volatile 
contaminants leaving a clean surface to receive the coating. 

Evaporating boats were fabricated from 5-mil tungsten ribbon and, when 
cleaned by preheating before use, proved to be satisfactory for the eva= 


poration of both gold and Si90 Temperature of the S10. evaporant was 


Do 
monitored with an optical pyrometer to maintain an indicated evaporating 


temperature of 2130°C; at this temperature 4.0 mg of $i0, evaporated in 
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about two minutes. The temperature is not corrected for departure of 


the evaporant from an ideal black body inasmuch as reproducibility and 


——_—_—_———_—sss 


temperature control were desired rather than an accurate measurement of 
true surface temperature. During the first five to ten seconds of all 
e€vaporations the evaporant was intercepted by a shutter, since the first | 
part of an evaporant to boil off is known to carry a higher percentage 
of any contaminants present. Following evaporation, cooling of the 
heater and targets was accelerated by admitting pure helium to the chamber 
through a liquid air trap. 

In practice, 2 gm or more of gold were evaporated at a distance of 
5 cm from the target backings and exposed parts of target assemblies to 


produce the desired plating thickness. 4.1 mg of SiO, were evaporated at 


2 
the same distance to give targets approximately 5-kev thick to 2-Mev pro- 
tons for use in the excitation function runs. For the targets used to 
obtain gamma-ray spectra, and for angular correlation runs at 1470 kev, 

5.2 mg of Sid, were evaporated to obtain targets 8-kev thick to 1470-kev 
protons. 

Considerable care was used in an attempt to reduce the contamination 
by oe. Evaporation boats were preheated to 2200°C before use, targets 
were stored in a vacuum (until this was found not to be of value), target 
chamber pressure was maintained at om Torr during runs, and possible 
sources of contamination in the vacuum systems such as Teflon were a- 
voided. In preliminary bombardments at the es resonance energy of 1373 kev, 
the gold-plated copper target backings, which were washed in either pure 
ethanol or menthanol and rinsed in distilled water just prior to use, were 
found to be freer from ce than backings of either Ta or W which had been 
cleaned by various methods. In spite of these precautions all targets 
which were prepared proved to be contaminated to some degree by fe 
2. Excitation function measurements 

The excitation function of the siine Goa Te oe reaction has been 
measured by recording single gamma-ray yield and gamma-gamma coincidence 
yield at accelerator energy settings over a range of proton energies from 
1420 kev to 2170 kev. The 2 in. x 2 in. NaI(Tl) scintillation detectors 
were used for this purpose. The output of one of the detectors was counted 
to obtain the gamma-ray yield at each proton energy, and the number of 
events in which gamma rays were detected simultaneously in each crystal was 


recorded by the use of the slow-fast coincidence system adjusted to a re- 


solving time of 25 ns. 
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The two detectors were positioned in geometry B at angles of 0° and 
180° relative to the vertical direction and at a distance of 0.6 in. from 
the beam spot. Integral discriminators in both channels were set so that 
counts would be recorded for gamma rays with energies exceeding 0.65 Mev, 
a setting chosen to exclude the 0.51-Mev annihilation radiation arising 
from the positron decay of Be while accepting the full energy peaks of 
the 0.68 and 0.71 Mev radiations which deexcite the first two excited 
states in Peo: The 7-microampere proton beam was collected by the cur- 
rent integrator which closed a beam shutter and stopped counters after 
collection of 84 microcoulombs at each energy step of approximately 3 kev. 
Background counts were made every 50 kev with the beam intercepted by a 
shutter in the beam tube outside of the shielded target area. 

An auxiliary experiment was conducted to determine the energy cali- 
bration for the excitation curves. In this experiment the frequency of 
the NMR signal was measured to within one part in Tom as the 992-kev 
aluminum resonance 2 was traversed using a freshly prepared target about 
10 kev in thickness. Runs with mass 2 beams were also made to check this 
calibration at a generator energy of 1984 kev. Immediately after cali- 
bration, the ee resonance energies above 1970 kev and the energy of 
the resonance at 1506 kev were carefully measured in a similar manner. 
Using these points as references it was then possible to calibrate the 
entire excitation curve in terms of magnet current, which was determined 
to a precision of one part in 10°? at each energy step. This procedure 
gave a better calibration than that obtained using the fluxmeter logging 
dial readings, since calibration curves for converting these readings to 
magnetic field could not be read with adequate precision. 

3. Measurement of saturation activity 

In order to definitely identify resonances due to the ce (pee 
reaction, measurements were made at the resonances above 19/70 kev to 
determine the half life of decay radiation following bombardment. In 
these measurements Sha targets were bombarded for ten minutes, then ree 
moved to a small cave of 2-in. lead brick where counting was started 
within 60 sec after cessation of bombardment. The positron annihilation 
radiation following beta decay of rae was detected by a 3 in. x 3 in. 
crystal. Total counts above the noise threshold were recorded by the 
pulse height analyzer operated in a gross counting mode. In this mode 


the channel in which all counts are recorded is automatically shifted 
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ahead after a preset counting time. After counting for 3.1 sec per 
channel for a 6.46-minute total counting time, the target was removed 
from the cave and a background run made with the analyzer in the sub- 
tract mode. 

Difficulty which was experienced in obtaining gamma-ray spectra 


ao 


recognizably due to Si at the upper resonances as a result of masking’ 
radiation from the oan O” reaction caused some concern that the 
peaks observed in the excitation curve might have been due to other cone 
taminants. Observance of the 2.5-min half life of pee following bombard- 
ment could not entirely exclude this possibility since it is conceivable 
for p?? to be formed in a nonresonant direct capture process. For this 
reason it was desirable to compare the activity of the target observed 
after bombardment at energies corresponding to the peaks in the yield 
curves to the activity observed after bombardment at off-peak energies. 
In an experiment to do this, a proton bombardment time of ten minutes to 
obtain 94 percent of saturation activity was followed by a ten-minute 
counting period during which the pulse height spectrum was obtained. 
The target was then allowed to decay for an additional ten minutes before 
bombardment was resumed at the next energy. The number of counts observed 
in the 0.51l-Mev annihilation radiation peak then served as a measure of 
the amount of yee formed at each energy. 
4. Gamma-ray spectra 

All gamma-ray spectra were obtained using 3 in. x 3 in. Nal(Tl) 
crystals with no special shielding nor collimation. For most spectra 
the pulse height analyzer usable range extended from about channel 7 
to channel 127 corresponding to gamma-ray energies from 0.7 to 7.7 Mev. 
For observation of annihilation radiation a higher gain was used so that 
these channels corresponded to energies from 0.18 to 1.92 Mev. The 
analysis of spectra for intensities was performed with a digital computer 
as described in Appendix 2. The particular method finally used to sub- 
tract background radiation from spectra varied according to the resonance 
being considered. When possible, backgrounds taken slightly above or be= 
low resonances were subtracted from spectra. This was not very successful 
for the resonances above 1950 kev where es contamination gave rise to the 
He (p.ome0-° reaction. This reaction is strongly resonant at 2030 kev 
with a width at half maximum of 120 ea which effectively blankets the 


higher proton energy range considered in the present work. The 6-7 Mev 


Zi 





(ts @e 
(ete 4 


gamma rays occurring in this reaction are present in relative proportions 
which vary considerably with bombarding energy so that off-resonance back- 
ground subtraction is not accurate. Bombardment of blank targets which 


had been processed identically to SiO, targets but using an empty eva- 


Z 
Poration boat gave a background measurement which was also rather unre- 
liable because of individual variations among targets. A method of sube- 
tracting the background measured with the accelerator off was also used, 
in which case background radiation of specific energies due to contaminants 
was subtracted out in the spectrum analysis. 

In order to circumvent radiation due to target contaminants efforts 
Were exerted to obtainsum-coincidence spectra in the energy region above 
1950 kev as a means of getting at decay schemes. The sum-coincidence 
spectra were, in general, measured with the slow-fast coincidence require=- 
ment also imposed, and with a coincidence resolving time of 110 ns. These 
runs were usually made in geometry B with counters positioned 1.6 in. from 
the beam spot at a 90° angle relative to each other. A 5/1l6~in. lead sheet 
between the counters served to minimize coincidences due to scattering or 
pair production in one crystal with the lost radiation being captured in the 
other crystal. The differential discriminator following the pulse adder was 
set with its lower level at the resonance energy and was adjusted to pass 
pulses in a band about 0.6 Mev in width. Settings were made by adjusting 
the discriminator so that the observed count rate of accepted pulses was 
equal to half the repetition rate of a test pulse. The proper test pulse 
amplitude was obtained by setting the pulser gain so that counts fell into 
the pulse height analyzer channel corresponding to the desired energy. 
settings of integral discriminators in the two singles channels were adjusted 
to 0.62 Mev in a similiar manner. The pulse height analyzer was calibrated 
with oe” and @°” sources and with the 5.3l-Mev gamma ray which dominates 
the setey Pe” spectrum at the 414-kev resonance. It was found that the 
Dumont 6363 photomultiplier tubes experienced marked gain shifts as a function 
of count rate amounting to several percent, an effect which has been attri- 
buted to changes in the dynode supnaceseee Attempts to allow this shift 
to stabilize before commencing counts were only partially effective. The 
overall result of this phenomenon was to vitiate the results of a large 


number of sumecoincidence spectra. 
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5. Measurements at 1470 kev 

Since the 4.50-Mev level in p20 has been shown to be most strongly 
excited at the 1470-kev resonance, this proton energy was a logical choice 
for angular correlation measurements which could lead to determination of 
the spineparity for the level. The resonance is rather weak so it was de- 
Cided to take five sets of data to improve the counting statistics and to 
Permit a check of the consistency of results. Each set of data consisted 
of five spectra taken with the movable detector set to angles of 0°, a0 
45°, 60°, and 90° relative to the proton beam. In geometry A with the 
3 in. x 3 in. movable detector positioned with its face 4.125 in. from 
the target, runs were made using a 7-microampere beam for counting times 
of approximately one hour at each angular setting. The same target was 
used for each set of rums and the order of the angular settings was varied 
in different sets so as to minimize differences in the gamma-ray spectra 
due to any contaminant buildup or target erosion. The fixed 3 in. x 3 in. 
detector was positioned at 270° relative to the beam and 4.125 in. from 
the target for use as a monitor counter. Each run was terminated when 
30,000 monitor counts had been recorded. A 2 in. x 2 in. detector was 
also used as a secondary monitor counter. The integral discriminator in 
each monitor channel was adjusted to a valley in the gamma-ray spectrum 
at 4.75 Mev so that the effect of gain shifts in these channels would be 
minimized. Before and after each run the discriminator levels were checked 
and adjusted with the aid of a ca source which was precisely and repro- 
ducibly located relative to the counters. For this check the gain in the 
monitor channels was increased by a factor of 4 by means of the gain mul- 
tiplier switch so that the discriminator setting was equivalent to 1.14 Mev. 
At this setting, near the 1.1/7°Mev peak in oe the count rate observed 
with the co source was quite sensitive to the discriminator setting so 
that accurate adjustment was possible. 

The unfolding of gamma-ray spectra and the determination of angular 
distribution coefficients was programed and analyzed on a CDC-1604 computer 
as described in Appendix 2. The radiations which have been used as spectro= 
meter calibration curves in this analysis are listed in Table l. These 
spectra were measured in experimental geometry A, backgrounds were sube 
tracted, and any remaining contaminating gamma rays were subtracted off. 
Each of these monoenergetic spectrograms was plotted, hand smoothed, and 


redivided into 525 channels to reduce any interpolation errors in the 
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TABLE 1 


Gamma Rays Used for Spectrometer Calibration 


E lev) Radiation Source SSD, 
0.662 Come : 

1.63 Rae (ene, one a2 
2.37 C “(outa 460 

2a 5 ll © 

oa Carne. 1698 
aaa BeLGpereiees 675 

ol Sis PS” ae 
6.135 Era Or. 1374 
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subsequent analysis. Energy calibration of each spectrum was made by 
equating distance in channels between the pair peak and the photopeak 

to 1.022 Mev. For the three lowest energy calibration curves, for which 
pair peaks are non-existent or poorly defined, an equivalent calibration 
was computed based on the pulse height analyzer calibration as determined 
in an auxiliary experiment. 

All spectra were read out of the pulse height analyzer memory photo- 
graphically. These data were then transcribed to punched cards and printed 
listings of the cards were carefully proofread against the photographs. 

Gain drifts did not prove troublesome at the lower count rates ob= 
served during correlation runs at 1470 kev. To compensate for differences 
in calibration of the pulse height analyzer between successive runs, which 
were of the order of one percent, each spectrum was self-calibrated by the 


peaks at 7.00 Mev and at 1.45 Mev. 
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IV EXPERIMENTAL RESULTS 


1. Excitation Function 

The excitation function for een as measured by gamma-ray 
yield is shown in Fig. 8, and as measured by gamma-gamma coincidence 
yield in Fig. 9. The energies at which resonances are observed are listed 
in Table 2, together with some of the results which have been obtained at 
proton energies below 1900 kev by other investigators. 

In arriving at the uncertainty of 0.20% kev for the resonance energies, 
consideration has been given to errors arising from instability and inhomo- 
geneity of the field of the analyzing magnet, the proton beam spread, pos-= 
sible surface contaminants, calibration drift, and the 0.5-kev estimated 
error for the 992-kev aluminum calibration point. 

Particular attention was given to the resonances above 1900 kev which 
have not previously been reported. The results of determinations of the 
half life of the post-bombardment activity of the targets at these energies 
are given in Table 3. The measurements are seen to agree with values for 
the half life of a reported in the literature, which range from 2.51 + 
.O2 min to 2.62 + .02 ing so that radiative capture by Sad is clearly 
demonstrated. Table 4 shows the intensity of positron Aneel Vecion 
radiation observed after saturation at selected bombardment energies. 

These intensity values show distinct maxima at 1975, 2033, 2075, and 2117 
kev and show that BP is formed in a resonant reaction rather than by a 
non-resonant process. The peaks in the gamma-ray yield curves at these 
energies can therefore be positively ascribed to resonant proton reaction 
with ane 

The yield of 0.511-Mev radiation after bombardment at energies between 
resonances is probably due to a combination of off-resonant reaction with 


see and to the 66-sec positron activity of ud following the non-resonant 


cms reaction. This hypothesis is supported by half life measure- 
ments of post=bombardment activity of the targets at two inter-=resonance 
Points at 1955 and 2007 kev. At those energies "half~lives'' of the order 
of 1.6 + .5 min are observed with the count rate 2 min after bombardment 


at 1955 kev a factor of ten lower than the count rate 2 min after bombard- 


ment at 1975 kev. 


26 





any 
- 
> 


















‘— oO — ie 
| ae 


bed 


i a 
oe 


































eh ne a 

BLL. ODA AL ml 

Co = —— es = 
—— — © Le a ee 
. | — —- a *"P & 
ee -_— oo ea 














——_ — °° = 
> am ~The 
7 a ont 











= &- => SS _ SS 
> 7 ~~... é 
a = 
a <> 
—— 
=> a> ‘Ga = SS ei a= 
————lC Ol kl Se ea e =a — 


ool2 


Q1TS1A 


0002 


AVY-VAWWYVS "8 SYuNndia 


(W3LSAS 8V1)49xy NI AOYSINS NOLOYd 
006! 008! oot 0091 


AW S90 (243 
ocd (2°49) 62! S 


0OSI 





oo! 


002 


oof 


00v 


GWOINOIOYXOIW Yad SLNNOD 


a7 





QTaIA JINIAGIONIOO 6 JuNdSI4 


(W3LSAS 9V71) 494 NI ASYSNZ NOLOYd 
oolz 0002 006! oos! 002! 009! oos! 


ocd (2*A) eS 


—a + 


ot 


Sl 


02 


Se 


BWOINOIJOMOIW Yad SLNNOD 


28 





TABLE 2 
, 29 30 , 
Resonances Observed in the Si “(p, %)P Reaction 


25 from 1450 to 2170. Values in kev. 


This work Cree Sear ondo lena Teyeoe | cepa = Remarks 
1470 1479 1479 1470 

1506 1515 1512 1505 

1637 1648 1643 1635 1638 

(1646) (1653) (1647) (650) 5a*° 
1667 1671 1663 eee 

1688 1692 1693 1680 1685 

(1726) weak 
752 ZZ 1752 1749 

1775 1777 1775 1772 

(1803) 1811 1808 

1860 1857 Lg32 

(1933) weak 
1975 

(2007) 517° 
2033 

2075 

2117 

(2159) not identified 
+0. 20%, +2 +2 +3 Errors 
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TABLE 3 


Half Life of Post-bombardment Decay of pro 


Bombardment Energy (kev) Half life (min) 
1975 2.42 + .09 
2033 2.45 + .18 
2075 2.51 + .06 
ZU] 2.47 + .15 


30 





TABLE 4 


Intensity of Post-bombardment Annihilation Radiation 


Bombarding neenuaLi on. Relative 
Energy (kev) Energy in P (Mev) Intensity 
1506 7.04 19.3 
1950 ~ 1.0 
1975 7.49 4.1 
1990 2 1.0 
2033 iS 7.9 
2060 = 1.8 
2075 Ooo 12.4 
2100 - Loe 
2117 7.63 8. 1 
2140 - 2a 


* 
Based on Q = 5.582, reference 24. 


el 








2. Gamma-ray Spectra at 1975, 2033, 2075 and 2117 k2v Resonances 

Gamma-ray spectra were obtained at the four resonances at proton 
energies of 1975, 2033, 2075 and 2117 kev. Masking radiation from the 
resonant reaction with the contaminant pee dominates the upper energy 
region of these spectra so that a detailed analysis of the relatively 
Weaker radiation from the $i-"(p, SP reaction was not possible. 
However, the gamma-ray energies and intensities which could be deter- 
mined with reasonable certainty from stripping and least squares analysis 
for intensities are given in Table 5. The 6=7 Mev gamma rays due to pl? 
which are omitted in Table 5, were accounted for by inclusion in the 
least squares intensity analysis, since other background subtraction 
techniques did not work well at these energies as has been discussed 
in Chapter III. Also included in Table 5 are energies of sum-coinci- 
dence spectral peaks which were close to energies observed in the 
singles spectra. The good spectral resolution and discrimination 
against background radiation which are obtainable with the sum-coinci- 
dence method proved to be of only nominal assistance in analyzing the 
gamma-ray spectra. This was due in part to the relative weakness of 
radiative decay of these resonance levels in comparison with inelastic 
Scattering, and in part due to photomultiplier gain drifts which blurred 
peaks in the sum-coincidence spectra. However, it is apparent from the 
large number of low intensity peaks observed in the sum-coincidence 
spectra at the higher-energy resonances that decay schemes from these 
levels are considerably more complex than decay from the 1752-kev re- 
sonance leva? where sum-coincidence spectra were recorded for the 
Purpose of comparison. 

The strong 1.28-Mev gamma ray observed at 1975, 2033, and 2075 kev 
is interpreted as deexcitation of the pe compound nucleus by proton 
reemission to the first excited state of aioe at 1.277 Mera The high 
intensity of this gamma ray is convincing evidence that this inelastic 
scattering process predominates over radiative capture and that the 
128=Mev radiation is primarily responsible for the upper-energy peaks 
in the gamma-ray yield curve. Rather curiously, however, this gamma 
ray does not appear in the spectrum at the 211/7-kev resonance. 

The resonance at 2075 kev lies near the Stee oieae resonance 


5 
reported at 2090 kev which has a half width of 16 kev. The gamma-ray 
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TABLE 5 


Gamma-ray Intensities at Resonance 
Energies of 1975, 2033, 2075, and 2117 kev 





E_ (kev) Singles Spectra Sum-Coincidence Intensity (quanta/«C) 
1975 4.52 4.58 
4.19 4.07 
| 3.48 
: ee 3.40 
3.00 2.99 3 
2.44 2.50 3 
1.28 - Ly 
2033 5.10 - 8.6 
26 - £3740 
20/2) 7.99 - 14 
4.70 4.85 60 
2.98 2.95 35 
2220 223 aD 
1.36 - 303 
1.28 - 1360 
2117 5.10 - Z 
23 = 3 


Error in singles spectra Ey estimated to be 1% 
Error in sum-coincidence Ex estimated to be 2% 


Error in intensities estimated to be 30% 
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TABLE 6 


Relative Gamma-ray Intensities at the 1470-kev Resonance 


* 
Ey (Mev) Transition Observed Decay Scheme’ Ref 23 
Intensity Intensity 





7.00 R~>0O 25.4 25.4 39.6 
6.29 R> 2 9.6 9.6 673 
4.50 19~ 0 14.6 14.3 13.2 
4.14 3 a= 0 10.7 1G; 7 CUE al 
3.98 R->9 10.5 10.8 13 
3.45 not accounted for 2.0 ~ 4.0 
Sa.7 R~> 11 3.0 i 2.0 
Oml7 ll~>l es: 2.0 
3.07 R= 12 19.1 17.6 De? 
3.02 9->0 1.2 1.6 
2.86 R>13 ioe 13e7 13.0 
2.69 133 276 228 - 
2.50 R19 38.8 21.6 20.0 
2.50 12-3 lj7a6 Ly 
2.34 9—~>1 ye 16.9 Isc 
1.45 19> 9 not analyzed’ 723 6.8 
1.45 30 AVS) ey 
a7 1 2—>0 bi 9.6 21.8 
. 68 1—>0 m 18.4 


4 An intensity of approximately 34 is indicated for the 1.45 Mev radiation 
by observations at 90° with .6"' crystal face to target distance. 


Relative intensities shown are also approximate absolute intensities 
in quanta per microcoulomb. 
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spectrum at this resonance is characterized by a 4./76-Mev radiation 

to the ground state of ee =e It seems probable, therefore, that 

the 4.70-Mev gamma ray seen in the spectrum at 2075 kev is due to the 
ee reaction. The observed intensity is in agreement with this ex- 
planation. Although the presence of 2.94-Mev and 2.24-Mev gamma rays 
in this spectrum suggest a cascade through level 8 at 2.94-Mev in ae 
the 1.45-Mev gamma ray known to occur in this maceaders is not seen 
and sum-coincidence work has not been able to confirm this possibility. 
It is possible, however, that the gamma rays due to git? reactions could 
be responsible for the anisotropy of the radiation from the 4.76-Mev 
level in pee which has been observed by hence: since anisotropy is 
not expected to result from the spin assignment to this level of 1/2+ 
based on proton scattering experiments. 

The remaining weaker radiations seen at the four higher resonances 
are most probably due to radiative capture by a but the rather large 
errors, due chiefly to masking contaminant radiation, do not permit more 
than tentative postulation of decay modes. Gamma rays observed at the 
1975-kev resonance suggest cascades through level 19 at 4.50 Mev and 
through level 14 at 4.18 Mev. No conclusions can be drawn for the 
2033-kev resonance. At the 2075-kev resonance the 7.59-Mev radiation 
is readily recognized as a ground state transition but the other gamma 
rays are not identified. Finally, at the 2117-kev resonance, the two 
radiations observed nicely fit a decay through level 5 at 2.54 Mev which 
is known to decay 95% to the ground yee. 

3. 1470-kev Resonance Decay and Intensities 

The gamma-ray spectrum observed at the 1470-kev resonance with a 
crystal face to target distance of 0.6 in. and at an observation angle 
of 90° is shown in Fig. 10. This spectrum has been least squares fitted 
by gamma rays with energies required by the decay scheme proposed in ref- 
erence 23. Reanalyses with alternate energy choices in each instance gave 
a poorer fit in reassurance of the correctness of the decay modes proposed. 
This decay scheme, with intensities as found in the present analysis, is 
shown in Fig. 11. (Heavy dots indicate cascades confirmed in reference 23 
by coincidence spectroscopy.) The branching ratio shown for the 4.50-Mev 
level is that given in reference 22. The assumed correctness of this 
ratio then requires branching ratios of 93:7 for level 9 and 80:20 for 


level 13 in order to fit the intensities of gamma rays which Have been 
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observed. These ratios are in reasonably good agreement with ratios of 
89:11 reported for level Es and 88:12 reported for level ears 
Intensities of gamma rays observed in the angular correlation runs 
are shown in Table 6 along with intensities observed in reference 23 at 
55°, The intensities found in the correlation runs have been averaged 
over 4n steradians by weighting the mean values for each observation angle 
by the sine of that angle. Because the 1.46-Mev background radiation due 
to ey was found to be quite strong relative to the 1.45-Mev 
gamma ray observed in the correlation runs, no attempt was made to analyze 
the gamma radiation of this or lower energy. Although the background 
radiation subtracts out, the statistics of the remaining counts do not 
allow extraction of significant information. The observed intensities 
are seen to agree reasonably well with those of reference 23 except for 
the 7.00-Mev radiation which is found to be considerably weaker in the 
present work. Since the experimental intensities are based upon reported 
values of intrinsic peak Semiereneiese: rather than an absolute spectro- 


meter calibration, this discrepancy is not considered significant. 


4. Angular Correlations at the 1470-kev Resonance 


A. General 

The correlation function coefficients which have been found in this 
experiment are given in Table 7. The values given are weighted means of 
the values found in each of the five sets of data taken. The chi-square 
values for each correlation coefficient have been calculated and show the 
different sets to be consistent, under the acceptance requirement that the 
cumulative chi-square probability be equal to or less than 0.9. Combining 
intensities from each data set, then analyzing for the coefficients gives 
values practically identical to the mean values quoted. This agreement 
to some degree serves as a check that computational errors have not crept 
in. Each experimentally observed spectrum has also been plotted together 
with the corresponding fitted spectrum, then visually examined to determine 
that a reasonable fit was obtained in each case. Fig. 12 shows a typical 
example of the fit obtained for the correlation runs. 

The intensity analysis was performed both subtracting background 
radiation analytically and subtracting off-resonance background, with 


essentially the same results being obtained in each case. The standard 
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deviations in the mean coefficient values have been calculated from 
n 
2 2 
1/6 = Zs L/é i 


and have been compared with the values given by 


i z a 
n(n-1) 2 a aa) 


which agree closely. 


The most obvious feature of the data is the anisotropy of the 
4.50-Mev gamma ray which is incompatible with a spin of zero for this 
level. It is also noteworthy that the 4.50-Mev gamma ray is the only 
radiation to show a statistically significant P, (cos @) term (notation 
is defined in Appendix 1). Fig. 13 shows the correlation coefficient 
values for this gamma ray obtained from each set of data. [In this 
figure M denotes mean values and C denotes values obtained by com- 
bining intensities from each data set before computing coefficients. 
Fig. 14 shows the fit to the Legendre series both with and without P,, 


terms. Although the existence of the Py, term is not incontrovertible 
it permits a markedly better fit to the data and gives a smaller value 
of the parameter € which serves as a guide to the appropriate number of 
terms which should appear in the Legendre series. 

In order to deduce spin-parity values it is necessary to examine 


the possible particle and radiation mixing ratios for this reaction. 


The notation used herein to refer to these mixing ratios is as follows: 


51 = Proton-partial-wave mixing parameter 
52 = The multipole mixing ratio of an intermediate 
unobserved gamma ray 


§3 = The multipole mixing ratio of an observed 
gamma ray 


4 = Channel-spin mixing parameter 


Appendix 1 defines these parameters more fully and describes how they 
enter into the expression for the theoretical correlation function. For 


the ehie eye and Si- (peer (intermediate unobserved) reactions, 
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the correlation coefficients Ay/Ap and A,/Ao are describable as functions 
of a maximum of three mixing ratios. Level curves, or contours, of con- 
stant A,/Ao values may be plotted against any two of the mixing ratios. 
Actually the arctangents of the mixing ratios are used as variables for 
convenience since they remain finite. If a third mixing ratio appears, 
changing its value changes these contour plots only to the extent of 
changing the A,/Ag values for all contours by a multiplicative constant. 
Contour plots for this analysis have been programmed and plotted as 
described in Appendix 2. 

In the following discussion it will be assumed that only El radiation 
or E2/Ml mixtures are likely to be observed for transitions from levels 
of several Mev. This simplifying assumption is based upon transition 
probabilities estimated from the independent particle model oe and is com- 
monly made in analysis of gamma radiation following proton capture in 
light and intermediate nuclei. Although suppression or inhibition of di- 
pole radiation is possible because of the isotopic spin selection rule, 
AT = + 1 for self conjugate ffeie:, there are an adequate number of 
lower levels of both T = O and T = 1 in Po so that the resonance level 
can decay by either El or E2/Ml radiation without violating this rule. 

In analyzing the experimental correlation coefficient values, tables 
from reference 23 have been quite helpful. This reference tabulates the 
maximum and minimum values which may be assumed by A,/Ao and A, /Ao for a 


variety of postulated spin sequences. 


B. Resonance level spin-parity 

The observed A,/Ao terms of about -0.40 for the 7.00, 6.29, and 
3.98 Mev gamma rays, which all proceed to known 1+ levels, are found to 
be most useful in determining the spin-parity of the resonance level. 
Among possible values through 4 for this level, it is found that only 
1+, 2+, and 2- can lead to the A,/Ao coefficient measured for the 7.00-Mev 
ground state transition. The same conclusion is supported by the A,/Ap 
values observed for the 6.29 and 3.98 Mev radiations. 

If one inspects the contour plot, Fig. 15, for the 1+-—)1+ transition 
it is seen that a d~s proton-angular-momentum mixing ratio of over 0.90 


is required to produce the observed A,/A value for the 7.00-Mev gamma 


0 
ray. Since the relative penetrability of d waves to s waves is about 


3 x on such a mixture is not a likely choice. The A, /Ao term observed 
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for the 4.50-Mev radiation supports this conculsion since A, /Ag terms 
do not result from resonance spins smaller than 2. Finally, the value 
of A,/Ap observed for the 2.34-=Mev radiation will be shown to be income- 
patible with a 1+ resonance level spin. This radiation proceeds to the 
known 0+ level 1, at 0.68 Mev, from the known 1+ level 9, at 3.02 Mev. 
Level 9 is fed more than 50% from the resonance level and, secondarily, 
by the intermediate 1.48-Mev radiation from level 19. The A,/Ag value 
observed for this primary 3.98-Mev feeder restricts the ranges of oi 
and g, for the R-»9-®1 transition. When these ranges are delineated 
on a contour plot for a 1+ -~—»1+ —>0+ transition, it is found that A,/Ap 
for the second radiation is restricted to a range between -Q0.01 and 
+0.24. By an Me naaenne of the formalism for intermediate unobserved 
cascades given in Appendix 1 it is possible to calculate the angular 
distribution coefficients for the 9-®1 radiation which result from the 
R->19-> 9-1 cascade. It is found that, for a 1+ resonance level, 


the minimum possible value of A,/A occurs for an assumed spin of 2 for 


level 19. This minimum value has — calculated to be -0.1/75. Thus 
a 1+ resonance level leads to a minimum possible value for A,/Ag of 
about -0.093. Comparing this with the observed value of -0.24 + .03 for 
the 2.34-Mev gamma ray, 1+ can be ruled out as a possible spin-parity 
of the resonance level. 

Having narrowed the possible resonance spin values to 2+ or 2-, 
the contour plots may be examined for the 2+-—~1+ and 2- -1+ transitions, 
applicable to the 7.00-Mev radiation to the ground state. In the first 
case, channel-spin mixing may occur while in the second case proton- 
partial-wave mixing is permitted. For the 2+—-»1+ case Fig. 16a shows 
the contour plot with the experimental value of A,/Ag for the 7.00-Mev 
gamma ray delineated thereon by cross-hatching. Fig. 16b is a similar 
plot for A, Ag: The darker shading in Fig. 16b indicates the superposi- 
tion of cross~hatching on both plots and defines the ranges and combi- 
nations of mixing parameters which can give rise to the observed coeffi- 
cient values. For $5 -20 (i.e. Ml radiation), §, is constrained to be 
less than tan 50°. This must not be taken too seriously, however, 
since the A, /Ao contours are widely spaced at the top and bottom boundaries 
of the permitted region, so that a very slight increase in the error for 


the experimental A /Ay value would allow So) to range up to 90°. The small 


i 
region at &, tan(-67°) and oF ~ tan 49° also fits the observed corre- 
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lation. The plots for the 2--»1+ case, shown in Figs. 17a and 17b, are 
quite similar to the 2+—™1+ case and show that there is no restriction 


on f-p proton-wave mixing for § -~0 (El radiation). The two small per- 


mitted regions at Sys sytan 67° = be ruled out since they require an 
M2/El mixture. It is seen that an unlikely but small M2/El mixture with 
> 3 7 0-026 or greater is required for this transition for the case of 
resonance formation by pure p waves. As will be shown, the 2- resonance 
level assignment is also less probable than 2+ based on mixing ratios 
which are required by experimentally observed coefficients for the 4.50- 
Mev radiation. 


As a result of the above considerations, a spin-parity of 2+ is proe- 


posed for the resonance level with 2- as a less likely possibility. 


C. R-—-y13->0 cascade 

There are only two cascades which allow cross checking of mixing 
parameters as a method of eliminating possible spin sequences. These 
are the cascade through level 19 and that through level 13 both of which 
levels have undetermined spins. The latter cascade will be considered 
first. The existence of an A, /Ao term which has been observed in the dis- 
tribution of the 4.14-Mev radiation from level 13 at the 1686-kev aoeeaanee = 
is evidence that level 13 must have a spin greater than or equal to 2. Of 
the likely possible spins for this level, 2+, 3+, and 4-, only 2+ and 2- 


can result in the negative value of A,/A observed for the 4.14-Mev gamma 


ray. An examination of the contour ioe for the four sequences 2+ -» 2+ 1+ 
for the R-»13-%0 cascade shows compatibility in all cases with the ob- 
served coefficients for the 2.86-Mev gamma ray and the 4.14-Mev gamma ray. 
For the sequence 2+->%2-—-~1+ the parameter 8, is constrained to be larger 
than tan 70° but the other possible sequences are compatible with a full 
range of values for either a or Sy, as the case may be. The only finding 
is that, for all cases, the R-®13 transition must be relatively pure di- 
pole radiation. The observed coefficient values for the weak 2.69-Mev gamma 


ray, which also is believed to deexcite level 13, are also compatible with 


all four possible spin sequences with reasonable choices for ¥, 


D. R-»19->0 cascade 
Since no requirements on parameters applicable to this cascade are 


imposed by transitions already considered it becomes necessary to consider 
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the R -»19 —0 sequence with no prior restrictions on the deltas involved. 
We find, furthermore, that no information is forthcoming from the R-?}l19, 
2.50-Mev gamma-ray correlation coefficient since this radiation is indis- 
tinguishable from the 12-3, 2.48-Mev gamma ray. 


A survey of maxima and minima of A,/A coefficients shows, if we 


0 
allow 2- as a possible resonance level spin, that the following sequences 
for the R-»19->0 transition are consistent with the experimental co- 


efficients for the 4.50-Mev radiation: 


2+ 1+ 1+ 2+ 1l- 1+ 
2- 14+ 1+ 2- le 14+ 
2+ 24+ 1+ 2- 2- 1+ 
2~- 2+ 1+ 2~ 3- 1+ 
2+ 34+ 1+ 24+ 3- 1+ 


It has been found in the analysis that a better fit (see Fig. 14) 
to the data is given by the series including a Py term. Then, provided 
this term is not spurious, which is not to be excluded as a possibility, 
the list can be reduced by eliminating sequences not giving rise to 
A, 1A terms. 


The following possibilities then remain: 


2+ 2+ 1+ 
2=- 2+ 1+ 
2+ 3+ 1+ 
2+ 3= 1+ 
2- 3+ 1+ 


The last two transitions shown require E3/M2 radiation so, adhering to 
the hypothesized El and E2/Ml radiative transitions as practical candi- 
dates, the spin of the 4.50-Mev level is limited to either 2+ or 3+. Of 
the three remaining possibilities, the 2- —»2+-%1+ sequence is perhaps 
least likely. Figs. 19a and 19b show that this sequence for R-19->0 
requires an f to p wave mixing parameter of over 0.10 with a large E2/Ml 
mixing ratio of nearly 12.0. For a 2+ =»3+—-~1+4+ sequence, Figs. 20a and 


20b show that an E2/Ml ratio of over 1.7 is required for the 4.50-Mev 
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radiation in order to fit the data. Because large E2/Ml mixing para- 
meters are not commonly seen in the lighter nuclei, this sequence, though 
more likely than 2~—-2+-—1+, does not appear as probable as a 2+ 32+ -1+ 
sequence, as will be seen. For the 24+-32+-m1+ sequence Figs. 18a and 
18b show that for an assumed channel-spin mixing ratio of zero, an E2/Ml 
mixing ratio of 0.27 will fit the observed coefficients for the 4.50-Mev 
gamma ray with a considerable latitude permitted for on Alternate 
choices for S$), will also permit acceptable choices of mixing parameters 
for this sequence though in general 8, and &, will be required to be 
larger. 

Summarizing, level 19 is found to have a spin of 2+ or 3+ under the 


assumptions which have been made, with 2+ being more probable. 


E. Comparison with other results 

Concurrently with the present work, an investigation of triple 
angular correlations at the 1470~kev resonance has been undertaken at 
Wright-Patterson Aeronautical wpewigenies. ™ The angular distribution 
of the 7.00-Mev gamma ray was also examined in this experiment. An A,/Ag 
value of -0.38+.02 was observed for this gamma ray in good agreement with 
the value shown in Table 7. Preliminary results of this experiment have 
indicated restriction of the resonance level spin to a value of 1 or 2 
and show that a zero spin for the 4.50-Mev level is not compatible with 
such a eeeenicenenw?- The preliminary analysis of the triple-correlation 
experiment is seen to be in agreement with the present work. It is to be 
hoped that the complete analysis of this work which is now in progress, 
will lead to a firmer indication of the spin-parity of the 4.50-Mev level 
Since the 2+ (3+) value proposed herein depends upon the somewhat debatable 
A, /Ao term observed and upon the hypothesis that only El or E2/Ml1 radiation 


is present. 


F. T= 1 character of level 19 

The prediction of 0+ for the spin-parity of level 19 is based on the 
premise that this level is the T = 1 analog of the 3.80-Mev level of si°’ 
in the isotopic spin multiplet Popmalaenaee A spin of O has been suggested 
for the 3.80-Mev level in ae based on observations of cascade gamma rays 


following the eecpeptysi-~ (tion. A spin of 0+ results also from 
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calculations made by Thankappen and Pandya of collective vibrational states 

in a weakly-coupled collective model of the AGcicwene The spin-parity of 
level 19 in pr? is thus expected both on experimental and theoretical grounds. 
The finding of non 0+ spin-parity for this level then appears to be good evi- 
dence that either the level is complex (e.g. a doublet) or that it is not the 
isotopic spin multiplet analog of the 4.80-Mev state in syed as has been ace 
cepted. Some evidence has been reported that the 3.80-Mev level in enue is 

a doublet, based upon observations of protons from the reactions | 

aa Gesi and Sconce: ae It is possible, therefore, that 
the analogous 4.50-Mev level in peo is a doublet, although no other evi- 
dence so suggests. The other possible explanation is that a nearby level 

in pe such as that at 4.42 Mev, may be the T = 1 level rather that the 
4.50 Mev level. This does not appear unlikely, since the T = 1 character 
attributed to the 4.50-Mev level is a probable rather than a firm assignment 
which was deduced from the investigation of the 532 (4a) pO reaction by Endt 
and pate” It is based upon a correction of the 3.80-Mev level in eo 

for the Coulomb energy and neutron-proton mass difference which yields a pre=- 
dicted energy of 4.47 Mev for the excitation energy in the Be nuclide. The 
assignment is strengthened by the fact that the alpha-particle group from 
s?*(d,a)P°o corresponding to a residual excitation in po of 4.501+ .010 Mev 
has only about half the intensity of neighboring groups. This has been inter- 
preted as being due to operation of the selection rule conserving isotopic 
Spin in this reaction. 

On the basis of the energy predicted for the T = 1 level in ae the 
4.42-Mev level is seen to be nearly as likely a candidate as the 4.50-Mev 
level, particularly since the energy prediction is not expected to be exact. 
Similarly in the T conservation argument, the observed lower intensity of 
the alpha group for the 4.50+-Mev level can hardly be considered conclusive 
evidence of the T = 1 character of the level, since some T = O levels are 
associated with alpha-particle groups of even lower intensity (e.g. the 
3.84-Mev level). Unfortunately, the 4.42-Mev level in Do is not known to 
be excited in radiative capture, so that spin-parity of this level will pro- 
bable have to be determined from a different reaction. 

5. Summary of Results 


Excitation curves for the jee reaction have been obtained over 


the proton energy region from 1420 through 2160 kev. Four previously unreported 
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resonances have been detected and confirmed at proton energies of 1975, 
2033, 2075, and 2117 kev with experimental uncertainties of +4 kev. 
Gamma-ray spectra at these resonances show that inelastic scattering pre= 
dominates over radiative capture, except at 2117 kev where no inelastic 
Scattering is evident. Some information concerning the radiative decay 
schemes at the 1975, 2075, and 2117 kev resonances was obtained by spectral 
analysis with the assistance of sum-coincidence spectra. 

The gamma-ray spectrum at the 1470-kev resonance has been analyzed 
and found to be in agreement with the decay scheme proposed by sigoney me 
The angular distributions of gamma rays observed at this resonance have 
been measured in a series of five experimental observations, each con- 
Sisting of gamma-ray spectra obtained at five angles relative to the pro- 
ton beam. Analysis of this data, under the assumption that El and E2/M1 
radiation deexcite me energy levels of several Mev, has shown that the 
1470-kev resonance level has a spin-parity of 2+ or, less probably, 2-. 
The spineparity of level 13 at 4.14 Mev has also been found to be limited 
to 2+ or 2-. Level 19 at 4.50 Mev was found not to have a spin-parity of 
O+ which had been eee ere Under the aboveementioned assumption, an 
apparent A, /A term observed in the angular distribution of this radiation 
leads to the conclusion that the spin-parity for level 19 is 2+ or, less 


probably, 3+. 
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VI APPENDICES 
Appendix 1 - Thecretical Description of Double Angular Correlations 


A. General 

There has been abundant literature coverage of the theory of angular 
correlations of successive nuclear radiations since Hamilton first pre- 
dicted the effect and developed theoretical expressions for successive 
emission of quanta. 40 A comprehensive review of the subject has been 
published by Biedenharn and Rose 46 and recent articles useful in exper- 
imental applications have been published by Gove 47 and by Biedenharm. 48 
The following notation is essentially that of Biedenharn. In all cases a 
primed symbol represents a larger value than the same symbol without a 
prime. 

In the case of a (p,) reaction in which the direction of a primary 
gamma ray deexciting a resonance level is observed relative to the direc- 
tion of an incident proton beam, the angular correlation is referred to 
as a double correlation. The function describing the relative probability 
of the gamma decay as a function of the angle between the beam direction 
and the emitted gamma ray can be written as an expression in the even 
Legendre polynomials: 

LD 


max 
w(e) = A (1)A, (2}B, (cose) » even (1) 
0 


ey = ' 
where A,(i) = a (LjL,)F, (L,L, 5,9) + 26a, (LL!) F, (L514 550) 


q 68 
+ §2a, (LIL, FB, (Luly 34d) 
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vay. i 
and F, (LL 3,0) = (yee 2 * aut 1) (+1) (20 +2) °® 
(LL'1 -1|vO)W(JSLL 32 5, ) 


The index h is related to the index i by: 


h = it-i¢l 


1 -1)v0) and W(JILL 52 J;) are respectively the 


Clebsch-Gordon coefficients and the Racah coefficients which are commonly 
used to describe systems involving coupling of angular momenta. The F 
functions, for which an extensive tabulation is available, 49 are defined 
in terms of a Clebsch-Gordon and a Racah coefficient by Biedenharn and 


Rose. 46 Arguments of the coefficients are angular momenta of the se- 


quence: ¥ 7" 
71 DB BL (2) 
Zz 2 
aa = channel spin = proton spin®@target spin 
L, Ly = interfering proton orbital angular momenta 
J = spin of resonance level 
L, BI = interfering gamma ray multipolarities 
Jo = spin of final state (not necessarily a ground state) 

S 1 = mixing parameter representing the ratio of the reduced 
matrix element for the transition involving Ly proton 
waves to that for Iy waves 

ie = mixing parameter representing the ratio of the amplitude 


of gamma rays of multipolarity Lo to those of multipo- 


larity L, 
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a_(LL') = cos(E L ~§ 5") » hare a's 1] 4 (3) 
is a particle parameter in which <a and 6, sare the phases at the nuclear 
surface of the particle partial wave functions of multipolarity L and lee 
The coherent interference between partial waves is expressed by the 
cos(@y “6') term. The factor g includes the coulomb phase shift due to 
electrostatic potential for charged particles and the phase shift occur- 
ring in "hard sphere" scattering. The a (LoL) terms for photons as 
opposed to particles are equal to unity. In the case where mixing of 
channel spins j, and j occur, W(@) may be expressed as an incoherent 
mixture: 

w(9) = Ww, (@)+ 67 Ws) (@) (4) 
where 5" is a mixing parameter expressing the ratio of probability of 
channel spin i entering the reaction to that of channel spin aa jy) 
The maximum order of Legendre polynomials appearing in (1), re is 
determined by conditions imposed on the W (abcd; ef) coefficients by con- 
servation of angular momentum which may be expressed by the so-called 
triad relationships. These require each of the triads of values (acf), 
(bdf), (abe), and (cde) to form a possible triangle. Restriction to the 
case of intermediate states of well defined parity (i.e. no interferring 
intermediate states of opposite parities) leads to the further simplifi- 
cation that Y may take on only even values. 

In the case of the reaction (p,x¥), in which x represents an unob- 
served intermediate gamma ray, the correlation is expressed as an inccher- 


ent mixture of multipoles occurring in the unobserved transition, i.e. in 
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the sequence: 


2\. 5 
4 fo J j (5) 
L 1 2 L 2 
@) Za 
lens eon oS elas ae) (6) 


where: 
Jo = channel spin 
Ly slo = interfering proton orbital angular momenta 
Jy = spin of resonance level 
L,,l, = interfering gamma ray multipolarities of unobserved 
radiation 


J. = spin of intermediate level 


2 
L,,L, - interfering gamma ray multipolarities of observed radiation 
J, = spin of final state 
&5 - ratio of intensity of the unobserved gamma ray of higher 


multipolarity to the intensity of the lower multipolarity 
gamma ray. 
The correlation function now contains an additional part expressing 
the connection between the two intermediate levels by the unobserved 


radiation as follows: 
DY max 


w(@) ys A (0)C, (1), (2)P, (cos9) (7) 


“td 


where the A, (i) are the same as in (1) and 
a uv By 
62) = (2) (20+ 1) (25g 41)) Fu(d, L,J53 3,5) 


It is seen that for the general double correlation case four mixing 
parameters may occur (barring the further complication of mixed inter- 


mediate states). 
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In the present investigation the Siz? target nucleus has a ground=- 
state spin of 1/2+, with the simplifying result that channel-spin 
mixing and proton-partial-wave mixing may not occur simultanously. The 
number of & 's is thus linited to a maximum of three. Explicitly, for 
resonance levels of 1+ ,2-, 34, 4-, etc. orbital angular inonentum mixing 
of two partial waves may occur but channel-spin mixing is prohibited. 
For the resonance level spins of opposite parities from those listed 
only channel=spin mixing is permitted. 

When three mixing parameters occur in the general case of inter-= 


mediate unobserved radiation (p,x8) with a spin sequence 


Lynd J 3 or j J J j (8) 
Orl 2 O 1 2 2 
Jo l,] “\L2 Loy '\ aimee 
the normalized correlation function can be expressed as 
w(6) = 1/Ap Dy AP, (cos®) (9) 


where each A [4 1s a product of the form: 


A/a = Hele ens.) t= (10} 
za (O at oyu yt ne 
with 4 = P5230 OF 25556 


and eres Bin gathe oie ea 


ZX 
ve (2+ 8%) 
(It is to be noted that the coefficient A,, appearing in (9) is a conven- 
tional representation of the more complete forms shown in (1) and (7). 
A distinction must also be made, of course, between the correlation funce 


tion W(®) and the Racah coefficient W(abcd; ef). 
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The coefficients a, b, and c can be expressed as follows: 


a, = aLo(lg +2) (Lodo) 
2L9(Lo + 1) -y@ +1) 
protcn , 
. 5, = sbolto + 1A) + 2))* aL 
. a Lo(Lo + 1)+ La(Lo +1) -ve+l) Fae{igto sory) my 
mixing 
. 20(Lo + 1) F (L! 
AL Agllg+ D = Av42) ee 
ayo 2 (LY (2d + 1) (239 + 1))]® WI EyJy5 JyJp) 
ungbs. 
gamma bo = O 
mixin 
Sg =e (2) [lam + 1) (2p + 1))® WH U5 y59450) 
a3 = FQ2J2d2) 
gamma | 
b3 = 2F(Lol2jaJq) 
mixing 
3 = F,, (LaJod9) 
es (La + 1) - 
es F, (Lg do71) 
at Qp(Lo + 1) = u(y + 1) Diane. 
chan. 
spin boy, ao 
mixing 
_, © 2Lo(lg+ 1) 8 
uh el eme ray wool 
Fy 5 cos( Gy = Sx) 
t2 = 43 > f, = 2 
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For the (p, ¥ ) case the only changes required are: 


eo O and 


pressed in terms of the coefficients Z and 2 


ficients’ 
al 
proton 
wave OL 
mixing 
Cy 


unobs, 
gamma 


mixing 


gamma 


mixing 


chane 
spin 
mixing 


| 


ay9 


v 


J, = J, = Jd. 


a 


The angular correlation function, as an alternative, may be ex~ 


00 


2 


(-1)90 Z(LoJ,Lody» Jor) 
2(=1)40 Z(LoJyLodys Jov) 


(-1)° 2(Egytosrs iow) 


at 
(-1) > W(J,Jpd,J03 Ly 2) 
0 


(1) 7H W( J, Ind, 403 1,2) 


(-1)"? 2, (Lodabodos Jo™) 


2(-1)"? 25 (hyd zLi'n5 3g%) 


(-1)°2 2, (LbJoboJns 5n%) 


(-1)"9 2(LodyLodys iov) 
0 


(-1) 50 ACTS REF 5gv) 
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ao i bia = 0, 


» instead of the F coef= 


(12) 





For (p,) reactions J, = J, = J, ayo = (2041), byp = 0, eyo £0 


il 
* (The two formalisms (11) and (12) are to some extent incompatible. 
Not only may the coherent mixing parameter phases differ by &1 20,04 


but also for Ly = QO the parameter &, doesn't represent the same matrix 
element ratio in both notations. The notation above, (12), has been used 
for most calculations in the present work. These expressions are in 


agreement with Sharp De 


except that the sign of b,,3 , for the (p,x%) 
case only, has been changed to agree with the part of Huby's correction 
applying to the observed gamma ray.209-4 For the "unnatural" radiation 
mixtures (ML/E2, E3/M2, etc.) tile sign of b,3 should be reversed. 

This sign convention permits the use of Moore's tabulated values <3 
where the same convention is used.) 

The Z and oa coefficients and Racah and Clebsch-Gordon coefficients 
have been tabulated by Sharp et al. = Using this tabulation Moore «2 
has prepared a convenient table of the a, b, and c coefficients applica-~ 
ble to radiative proton capture by gi°?, 

In order to fit experimental data on angular correlation coeffi- 
cients to theory, it is desirable to first eliminate as many spin se= 
quences as possible on the basis of any evidence which may be available. 
In cases where only two mixing parameters enter, contour plots 47 showe= 
ing contours of constant theoretical values of the coefficients A/Ao 

can be plotted against the inverse tangents of the two mixing parameters 
for all possible spin sequences. If three mixing parameters enter, one 
of them may be held constant for the purpose of plotting the correlation 
coefficient contours against the inverse tangents of the other two param=~ 
eters. It is seen from (10) that the contour values for different 


values of the fixed parameter differ only by a miltiplicative constant. 
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When experimental values of correlation coefficients are delineated on 
the contour plot, some spin sequences can be eliminated and the possible 
combinations of mixing parameters are restricted for other sequences. 
The observed values of several correlation coefficients in one decay 
scheme may then fix or highly restrict values of the spin-=-parities of 
resonance and intermediate levels, and determine the values of the mix-= 
ing parameters. 
B. Application to 1470-kev resonance 

For the resonance occurring at a laboratory proton energy of 1470 
kev, the phase factor cos (6, -§ 04) has been calculated using the 
expression Gr, =6,+ ?, where 6. represents the phase shift due to the 


point coulomb potential and py is the “hard sphere" phase shift. 


6,, is defined 42 by exp(2i€,) = PM +1+ in) (13) 
T(L+1-+im) 
' 2 
where ” = = which leads tos: 
Vv 


6i4276L = tan” (2) + tan” (74+) (14) 
The value for Pr defined by 
~ OC mskR) = iF, (m,kR) 
Gy (mq skR) + LF; (™ okR) 
where F and G are the regular and irregular solutions to the radial wave 


exp(2i $,) (15) 


equation for the coulomb potential. The value for p, can be found from 
the graphs of Sharp et al. 23 (these graphs actually list~@ ) and em-= 
ploying R = 1.54/36, as an estimate of the nuclear radius. These 


graphs have also been used to calculate the relative pentration factors 


Pe + G 
v,/Vp 2 = 11-2 (16) 
F + Gz 
which are useful as a guide to the expected mixing parameter of inter= 


fering partial proton waves. Results of these calculations appear in 


Tables 8 and 9. 


71 





TABLE 8 


Relative Penetration Factors for 14/70-Kev Protons 


d/s f/p g/d 
ee oO eee i.3°x 10” 
TABLE 9 


Coulomb Phase Factors for 1470°-Kev Protons 


€.-§, €,-€, €,-€, 


O O Oo 


104 73 56 


ee 





Appendix 2 - Computer Programs and Analysis 


A. General 

In order to analyze experimental data in the form of gamma-ray 
spectra recorded by the 128-channel analyzer, two computer programs 
were developed. The first and more extensive program entitled program 
SPECFIT computes experimental correlation coefficients, while program 
CONTOUR is designed to utilize a California Computer Products Model 
565 graph plotter as an output medium to automatically plot contours 
of constant angular correlation coefficient values against the inverse 
tangents of two mixing parameters. Both programs are written in fortran 


language and were run on a CDC-1604 digital computer. 


B. Program SPECFIT 

The output pulses of a photomultiplier tube viewing a scintillation 
crystal upon which monoenergetic photons impinge do not all have the same 
voltage. Instead, the three fundamental processes of photoelectric ab- 
sorption, Compton scattering and pair production, which occur in the 
crystal with probabilities dependent upon gamma-eray energy, lead to a 
response pulse-spectrum which is characteristic of each particular gamma= 
ray energy. Experimental geometry, crystal size, shielding, various line- 
broadening effects, and characteristics of individual crystal-photomulti- 
plier tube combinations influence the response also, so that response 
functions are in general unique for each spectrometer in a given experi- 
mental arrangement. 

An experimentally observed spectrum is, of course, usually composed 
of a number of.gamma rays, the spectrum being the sum of the constituent 
radiations each multiplied by its intensity. The problem of deciphering 
such an observed spectrum to obtain the energies and intensities of com- 
ponent gamma rays is of fundamental importance in gamma ray spectroscopy 
and has not yet been satisfactorily solved in complete generality. 

Using notation, which is appropriate for the discrete pulse height 
intervals into which multi-channel analyzers sort pulses, and also ap- 
propriate for the discrete number of gamma ray energies composing re= 


sonance deexcitation spectra, one can write; 


Ss, = 2, RI, 


™m 
j=l 


us 





where S. = observed number of counts in channel i 


a 
R. = probability that a photon of energy E. 
J originating at a particular source 
will produce a pulse falling in pulse 
height channel i (usually referred to 
as the response matrix) 
o = the number of photons of energy Bs 


originating at the source 


m = the number of gamma ray energies which 
occur in the spectrum under consideration 


A rather fundamental and obvious approach to this problem is to 
assume no knowledge of gamma-ray energies is available, divide the gamma 
Yay energy range into a number of discrete intervals equal to the number 


of channels, then solve for: 


In fact, the first attempt to unfold experimental spectra in the 
present work utilized a modification of a apestaaee which employed an 
iterative technique developed by Se oereline to solve the above equation. 
This attempt failed for two reasons which are rather instructive. First, 
the analytical expression used to approximate the response matrix did not 
adequately represent the true spectrometer response, and second, the 
technique proved inadequate with experimental counting statistics attain- 
able to handle the analysis of approximately 100 channels. 

In order to avoid the first problem resort was had to the technique 
of interpolating between experimentally-determined moncenergetic cali- 
bration curves in order to obtain the response matrix. In order to over= 
come the second objection it was decided to use information on the gamma 
ray energies available from preliminary "'stripping"’ of spectra and from 
other sources. Having thus reduced the number of gamma rays to be con= 
Sidered, a linear combination of the selected response matrix elements 
was then fitted to the observed spectrum by a least squares process. 


This procedure was followed by another least squares analysis to obtain 
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correlation function coefficients from the gamma ray intensities cal- 
culated for a series of runs at different experimental angles. The 
steps in the computation procedure in program SPECFIT as finally de- 
veloped are as follows: 

1. Load data. (See section D of this Appendix) 

2. Subtract any constant background from calibration spectra. This 
background represents the Compton "tails" of higher energy gamma rays 
contaminating certain of the calibration gamma rays. 

3. Normalize each calibration spectrum to a photopeak height of one 
unit. 

4. Include a straight line in the response matrix if this option is 
selected. This step is essentially a refinement of step 2. 

5. Correct the input and background spectra for dead time of the pulse 
height analyzer. For this purpose the dead time given by the manufacturer 
is 38 + 0.5i microseconds where i is the pulse height channel number. 


This leads to an expression for corrected counts: 





N. 
| es 
= 127 


1 - 1/T >, N, (38 + 0.51)*10°° 
2 


the number of counts observed in channel i 


where N. 
i 


T = counting time in seconds 


6. Calculate energies in Mev corresponding to each channel number of 
the input spectrum. 
7. Adjust energy scale and number of channels of background spectrum 
to match the input spectrum. Multiply by ratio of run times to get the 
proper background to be subtracted from the input spectrum and perform 
the subtraction. 
8. Using the monoenergetic calibration spectra, calculate a matrix of 
response curves for the gamma rays which the spectrum is to be decomposed 
into. This calculation is performed in a subroutine entitled RESPONS as 
follows: 

a. Find the two calibration curves bracketing a gamma ray energy 


desired as a response matrix element. 
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Shift the lower-energy calibration curve found in a. so that 

its photopeak falls in the same channel as the higher energy 
photopeak. Fill in bottom of shifted curve as determined by 
the slope between the bottom two unshifted channels. Match 
energy calibration of the lower-energy curve to the upper. 

(For this step, and elsewhere where interpolation is necessary, 
cubic interpolation is performed using the four nearest points 
to the desired value.) 

Interpolate linearly between each pair of points of the 
bracketing calibration curves to obtain a response curve for 

the desired gamma-ray energy. (A "curve" in this description 
is meant to imply the vector consisting of a 1 X n array of 
values, n being the number of channels used in the description). 
Using the input spectrum channel energies calculated in 6., 
interpolate to get a response curve with channels and channel 
energies corresponding to those of the input data. 

Normalize the response curve for the given gamma-ray energy so 
that the ordinate corresponding to each channel represents the 
probability that one photon leaving the source will produce a 
count in that channel. For this calculation the experimentally- 
determined intrinsic peak efficiency data of Lape sheen een 
used to derive the following approximate expression for intrinsic 


peak efficiency as a function of source distance: 


I 


€ = ar * 


-165 + .Q067h 


> 
l 


where 


©9699 - .00918h 


” 
It 


h = source to crystal distance in centimeters 


E = gamma ray energy in Mev 


Assuming an accurately gaussian photopeak, and that the crystal 
resolution is approximately inversely proportional to the square 
root of gamma ray energy, (an assumption which was experimentally 


verified) the ratio of photopeak height to area can be expressed: 


1 2 | Gn2 1 
C= Ee ee = 
6Y 2x WR WE» 2020 .866WR YE 
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fe 


where R is the usual resolution definition of peak width at 
one half maximum amplitude divided by energy for the come 
gamma ray, and W is the spectrometer calibration in pulse 
height channels per Mev. When @ is multiplied by the appro- 
priate solid angle fraction and attenuation coefficient, the 
value obtained represents the probability of recording one 
count in the photopeak channel when one photon originates at 


the source. 


Return to main routine. 


9. Fit the experimental spectrum with a linear combination of the selected 


monoenergetic response curves, the coefficients representing absolute in- 


tensities because of the method of normalizing the response matrix. The 


weighted sum of the squares of the differnces between observed channel 


counts and fitted counts is minimized in order to arrive at this fit. 


The normal equations resulting from the minimization are: 


S.R..W. = R..R..I.w.=0 ike } 
2 eal cos) ya Jee ae <3 
j iy 


where oe = Observed counts in channel j 


Rap = Response matrix element giving 
J probability of a count in channel 
j due to one photon of energy E 
aa , k 
originating at the source 


w. = Weight assigned to the squared 
J error in channel j 
I. = Number of photons of energy E, 


which leave the source during 
the counting period 


j summations are over the channels for 
which a fit is desired 


i summation is over the number of gamma 
ray energies E. used to fit the ob- 


served spectrum 


= 


For calculation the 1. are expressed: I. - 2, Avs By 
with Ala = 2 RR KY; and By. = ae 4 
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A standard subroutine (GAUSS3) is used to invert A. Mention should be 
made of the weights e which have been used. Normally, for a Poisson 
distribution the weight is taken as Ke where the number of observed 
counts in a channel are used as an estimate of the variance 67. 

ae: tener who has pointed out that this procedure overweights statisti- 
cally low counts, identifies the Poisson distribution variance, which 
rightly is the mean number of counts, with the fitted number of counts 
for each channel. This procedure is adopted here with a modification to 
account for statistical error arising from the subtraction of background 
counts. We have Ss. = S, - b,, where S, are observed channel counts cor- 


rected for dead time and b. are background counts to be subtracted. Then, 
ae 2 2 
6°(S!) =6°(S,) + 6 (b,) 
-, 4 
7S, + (T/T,)b, =S. + (1 + T/T, )b, 


where T = counting time for data and Ty = counting time for background 
run. Identifying Si with the computed number of counts, Sif: then gives 
Si + (1 + T/T, )b, for the channel variances. To avoid overweighting of 
very small values of fitted counts an arbitrary minimum value of 6 is 
chosen to be two per cent of the average value of 6 over the fitted 
channels. Since no fitted count values are initially available, an ini- 


tial weight of one is assigned to each channel and the solution is then 


iterated to convergence, stopping when 


2 
Ss; (81) is smaller than 107° 


i 
(This series of iterations, usually five in the present work, accounts 
for most of the program running time and can be abbreviated for prelimi- 
nary analyses; see reference 57 in this regard.) é1,/1, is the fractional 
change in the computed intensity of the ia gamma ray between successive 
iterations. 
10. The computed intensities are examined for algebraic sign. Since 
only positive intensities have a physical basis any negative intensities 
are discarded and the problem resolved to this point. 
11. For each gamma ray calculate the fitted total number of counts in the 


channels above the channel representing 40 per cent of the photopeak energy. 
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(This is an arbitrary calculation used only for checking results obtained 
in manual unfolding). The fractional error assigned to these counts is 

set equal to the fractional error in the corresponding intensity. The 
statistical variance estimators of the intensities have been taken to be 
the diagonal elements of the inverse of the matrix A in a, In addition 
to statistical intensity errors there is an error due to gain and discrimi- 
nator drifts in the monitor counts channel which must be accounted for as 
well as error due to eccentricity of the beam spot on the target. These 
fractional errors have been estimated to be 3.3 per cent and 1.2 per cent 
respectively in the present experiment. The three errors are combined 


according to?; 


1h >, /6 , 


J 


Two other possible sources of error should be noted here. First, the 
errors in the experimentally determined response matrix, which can con- 
siderably complicate the error calculations for the intensities, are kept 
to a minimum by obtaining good counting statistics in the calibration 
spectra. The other source of error is gain drift which tends to smear 
spectral peaks. Fortunately, the main effect of such drift seems to be 

an increase in the mean squared errors in the curve fitting process and 
intensity determinations are apparently not too sensitive to the effect. 
12. The intensities and their errors calculated in 11 are read out and 
stored for use in the subsequent determination of correlation coefficients 
and mean values. 

13. The absolute values of differences between observed and fitted counts 
are summed over fitted channels and divided by total input spectrum counts 
over the same range to give a measure of the fit which has been achieved 
in the analysis. This measure is then read out. 

14. The input and fitted spectra are plotted in a graph output. 

15. The above process is repeated for each run in the data set until the 
data for all the experimental angles has been processed. 

16. When the preceding analysis is completed for a set of data a sub- 
routine entitled LEASTSQ may be called, if appropriate, to compute the 
correlation coefficients. For each gamma ray in the spectrum a least 
squares fit of the counts computed at the experimental angles is made to 


the usual series: 


1 
W(@) = 23 AP, (cos Q@), i even 
i 
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where P. are the even Legendre polynomials. The fit is first made to 


Ay + AP, + A,P) then Ap + ADP, - 


to the series terminating with P 


(Provision has been made for fitting 
6 by simply changing the initial value 


of NP to 4; however, see reference 61 concerning accuracy of P, terms). 


6 
The computation, which follows the procedure given by Rose = proceeds 
as follows: 

a. Table gamma ray absorption ee aehenee in copper and cor- 
responding energies in Mev. Correct intensities and sigmas 
for attenuation in copper target backing. 

b. Set up matrix ah of Legendre polynomials for experimental 
angles. 

c. Do least squares calculation to find coefficients using the 
reciprocal of the estimated intensity variances from 11 as 
weights. 

d. Normalize coefficients and correct for the effect of the 
finite solid angle subtended by the scintillation crystal. 
These latter corrections are precomputed in an auxiliary 
program using the expressions given by row Estimated 


variances of the coefficient ratios, A,/A etc., are com- 


0? 
puted from the expression given by Reich et le (After 
computation variances were increased by a factor 1/6 é in 
all cases where the computed value of e* was less than unity.) 

e. Calculate anisotropy. 

f. Compute the cpacemignt This quantity, defined by Rose's 
analysis, is a measure of non-statistical errors in the data. 
If such errors are eliminated the value Bnee serves to in- 
dicate whether the appropriate number of polynomials has been 
used in the correlation function expression, with values of 
Saterescer than unity indicating an inappropriate choice. 

It is pertinent here that non-statistical errors have already 
been introduced into the problem in 11 above so that this 
criterion is valid only to the extent that an accurate esti- 
mate of the true variance of the counts has been achieved. 
The computed values of a in general are near unity and seem 
to verify that the estimated variances can be considered to 


be reasonably good. 
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g. Selected values of the normalized calculated correlation 
function, and values of a comparison theoretical function 
if desired, are read out, and similarly normalized experi- 
mental points are read out. (Experimental values are un- 
corrected for solid angle, the correction being applied to 
the theoretical function for comparison with experimental 
results.) 

h. Iterate for each gamma-ray energy. 

i. Reduce number of polynomials if orders higher than P, have 
been used and iterate solution. (The program also calculates 


coefficients for the cosine series: 


W(@) = 1/4) 5S Ay, [cos (@)]*" ) 


0 


~ Ms 


j. Return to main routine. 
16. Repeat entire analysis for as many sets of data as indicated. 
17. Combine coefficients to get weighted mean values and Chi square 
values if more than one set of data is analyzed. Read out these values. 
18. Combine intensities at each angle from all sets of data and re=- 


analyze for coefficients using this combined data set. 


C. Glossary of names of selected variables and constants appearing 


in program SPECFIT. An asterisk indicates an input quantity. 


ACI, J) The matric product of the response matrix with its 
transpose, weighted with the WT values. 


AINV(I, J) The inverse of A(I,J). 

ANISO(CT) Anisotropy for Legendre series solution (I = 1), 
and for cosine series solution (I = 2). 

ARGP(T) Vector containing cosines of the angles B3(TI). 

*B(L) A vector representing the product of the inverse 


of the response matrix and the observed spectrum. 


*BACK(TI) A number representing background counts to be 
subtracted from calibration curve I. 


BACKGND (1) The calculated background spectrum to be subtracted 
from an input spectrum. 


*BCHAN The lower of the channel pair described under TCHAN 
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FIN(1,J,K) 


FMU(I) 


*HCM 
HOLD (1) 


LPKS 


*MAXCHAN 


*MESH 


*MINT 


*MINTB 


MPKS 


MPKS 2 


*NBAOPT 


*NCUT 


NERFL 


*NHIEST 


*NLOWEST 


*NOCALB 


*NPHOTO 


Storage for FINIENS(1), data set J, run K. 


Mass absorption coefficient in copper for EMU(1) 
Mev gamma ray Units in cm4/gm. 


Distance from source to crystal face in centimeters. 
A temporary storage vector. 


Same as NPKS but not including any straight line 
component. (See NBAOPT) 


An upper channel limit in loading CALB data. Between 
MAXCHAN and MINT, CALB(I) values are filled in with 
the value of CALB(MAXCHAN). 

The number of channels in input spectrum. 

Number of channels in calibration spectra. 

Number of background channels. 

Initial value of NPKS for restoring the value. 


Initial value of LPKS for restoring the value. 


Background option selector. 


NBAOPT = O No background read in nor subtracted 
from SPEC 

NBAOPT = 1 Background spectrum is read in and 
SPEC is corrected for this background 

NBAOPT = 2 Same as NBAOPT = 1 and in addition a 


straight line is included in the 
response matrix. 


A channel number below which the graph output 
values are divided by a scale factor SCALE. 


An output of matrix inversion subroutine. A value 
of 1 indicates matrix is non singular, and 2 indicates 


a singular matrix. 


The upper channel of the range over which a fit is 
made in SPECFIT. 


The lower channel of the range over which a fit is 
made in SPECFIT. 


Number of calibration spectra. 


The photograph number of the observed spectrum; for 
identification only. 
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*NPKS 


*NRUNS 


NP 


*NPRINT 


*NSETS 


OMEGA (I,J) 


*PKSEP(1) 


* PHOT (I) 


*PKEN(I1) 


POLY (I,J) 


R(I, J) 


RED 


*RES 


SAN(1) 


SAVE(I,J) 


The number of spectral components with which the 
observed spectrum is to fitted. Program will increase 
this value by one if option is chosen to include a 
straight line in R (See NBAOPT) . 


Number of runs in a set of data. 


Number of Legendre polynomials to be fitted in wirse 
pass through LEASTSQ. 


Print suppression index 


NPRINT = 0 All printouts and graphs are made 

NPRINT = 1 Suppress printouts of spectrum values and 
cosine series coefficients 

NPRINT = 2 Same as 1 and also suppress graph outputs 

NPRINT = 3. Same as 2 with briefer more compact arrange~ 
ment of correlation coefficients 

NPRINT = 4 Same as 3 but graphs plotted 


Number of sets of data. 


Used for temporary storage in SPECFIT and diagonal 
elements used for weights in LEASTSQ. 


The number of channels corresponding to 1.022 Mev 
for energy calibration of CALB curves. Taken as 
channel distance from photo peak to pair peak for 
higher energy calibration gammas. 

The channel of calibration curve I in which the 
photopeak occurs. (May be an interpolation between 
two integral channels numbers). 


Energy in Mev of a spectral component. 


The matrix of Legendre polynomials of order 2(J-1) 
evaluated at angle B3(I). 


The response matrix. 

A normalizing factor for 3" x 3" unshielded Nal(TlL) 
crystal representing the probability that one photon 
at the source will produce one count in the photopeak 


channel. 


The proton energy in kev at the resonance being 
investigated (for identification only). 


Estimate of standard deviation of ANISO(I). 


An auxiliary matrix. 
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SBB(1) 


*SCALE 
SCOEFT(I) 
Sigems (1) 


SIGTENS (1) 


SIN(1I,J,K,) 


*S PEC (TI) 


STC(1,J,K) 


STORE(I,J) 


STSC(I,J,K) 
STW(1I,J) 
STSW(I, J) 
SUM7 (I) 
*TBACK 


TC(I) 


*TCHAN 


TEMP (1) 
*TOPEN 
TS (1) 
*TX60SEC 


WT(I) 


XI(I) 


Coefficient for renormalizing input spectrum with 
counts from an auxiliary monitor (optional; if not 
entered values are unity). 

A scale factor for graph output. See NCUT. 

The standard deviation in COEFT(I1). 

The calculated error in COUNTS(I). 


The error for calculated spectral intensity 
INTENS (I). 


Storage for errors in FIN(I1,J,K) 


The vector representing the counts in the observed 
spectrum. 


Storage for accumulating calculated correlation 
coefficients of order I, for gamma ray J and 


data set K. 


A temporary storage matrix for A while it is being 
inverted. 


Standard deviation of STC(1,J,K). 

Temporary storage. 

Temporary storage. 

Temporary storage. 

Duration of background run in units of seconds x 60 


The weighted mean value of correlation coefficient 
of order I. 


The upper of a pair of input spectrum channels to be 
used for sensing energy calibration of the input spectrum. 


Temporary storage. 

Energy of TCHAN in Mev. 

Standard deviation of TC(1) 

Duration of a run in units of seconds x 60. 


The weight assigned to channel I in least squares 
solution for intensities. 


The matric product of POLY and Bl weighted with the 
diagonal elements of OMEGA. 
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XS (1) A vector representing the energy in Mev of each 
channel of the input spectrum. 


Z(1) A temporary storage vector. 


*ZERO The number which if any diagonal element of matrix 
A is smaller than, A will be consideral singular. 
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D. DATA INPUT SEQUENCE FOR PROGRAM SPECFIT 


Number Number of 
Input of sets cards in set Input name Remarks 

1 ih 1 NSETS, NPRINT 4 

2 i 1 EL, ECC, HCM 

3 1 l NOCALB, MINT 

4  NOCALB 1 BACK(J), PHOT(J), ENCALB(J), 
MAXCHAN, PKSEP(J) 

MAXCHAN/7 CALB(I, J) 

5 1 1 NPKS, NBAOPT 

6 1 NPKS PKEN(1) 

qj 1 1 BSL, BEN, BPHOT, MINTB, l 
TBACK 

8 i MINTB/7 BIGT(1) 

9 1 it NRUNS, CTSNORA, NLOWEST, 
NHIEST, NCUT, SCALE, ZERO 2 

10 i 1 RES, NPHOTO 3 

11 1 1 MESH, TCHAN, TOPEN, BCHAN, 
BOTEN, B3(NRUNDEX), T60XSEC 

12 1 1 CTSNORB 

13 l MESH/7 SPEC (1) 

14 1 4 Graph data cards 5 

REMARKS : 


1.) Read in only for NBAOPT greater than l. 
2.) Input resumes here for each successive data set. 
3.) Input resumes here for each successive run within a data set. 


4.) All loading formats are multiples of F10.3 and I10 except that ZERO 
is El10.1 and TBACK is F20.3. 


5.) Blank cards may be loaded if graph option not used. 
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E. Description of program CONTOUR 


As described in appendix 1 the correlation function coefficients for 
(p,f ) and (p,x¥) intermediate unobserved reactions for spin 1/2 target 
nuclides can be described as a product of, at most, three functions of the 
form: a 

a aint Pin Oe t in Ok 
2 1 +60 


where i denotes the order of the coefficient and k is the index of 
the mixing parameter & . For computation the substitution S = tan Y 
is commonly used, and Weade”™ to ¢ =Q + 6 sin (+7) 


a+c 
Z 


yee 
B=+5 Eeeacentor| 





with oO 


x ca Gucee ) 


b/ |b] 


WM 
il 


where indices have been dropped for clarity. 
This program computes sets of possible values of two mixing para~ 
meters for various fixed values of correlation coefficients. Using a 
graph plotter output medium available through a locally developed sub- 
routine, these sets of values are plotted to form a set of contours 
along which the value of the correlation coefficient is constant. 
Steps in the computation are as follows: 
1. Read in data for one possible spin sequence to be examined. 
Compute values of @, B and & for the lowest order coefficient 
which occurs. If three mixing parameters occur simultaneously, 
read out the expression for the fixed g as a function of the 
fixed mixing parameter. (All loading formats are standardized 
to be F10.3 and 110 for floating and fixed point variables 
respectively.) 
2. Find and read out the maximum and minimum possible values of 


the coefficient and the corresponding values of the Y's. 
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Fix the coefficient at an initial value within its range. 
Assign values to one of the Y's in equal increments through- 
out the range -90° to +90° and solve for the corresponding 
values of the remaining Y, which in general will be double 
valued over the region selected. 

Read out pairs of Y values for the selected coefficient 
value both to the print medium and to the graph plotter. 
Repeat steps three and four for successively larger values 
of the coefficient until a specified number of contours has 
been plotted. 

Repeat steps two through five until contour plots for co= 
efficients of all orders permitted by the given spin sequence 


have been produced. 


F. Input variables for CONTOUR 


NSETS 


JF IND 


JSTEP 


JF IX 


NCONTS 


NPTS 


DELTA (JF IX) 


Number of spin sequences for which plots are to be 
made. 


Index of the mixing parameter to be solved for. (Y 
values on graph). 


l = proton partial wave mixing 


It 


unobserved gamma multipolarity mixing 


H} 


2 
3 observed gamma multipolarity mixing 
a 


channel spin mixing 


Index of mixing parameter assigned successive values 
(x values on graph) See JFIND. 


Index of mixing parameter to be held constant when 
three parameters occur simultaneously. (Set to five 


for (p, ¥) reaction where only two parameters occur.) 


Index specifying the highest order coefficient occurring 
in the spin sequence. (Order = 2¢eIMAX) 


Largest value which the index of mixing ratio assumes 
for use in reading input data. (Usually set at four.) 


Number of contours to be plotted on each graph. 


Number of pairs of points to be computed for each branch 
of one contour. 


Fixed value assigned to the third mixing parameter during 
a calculation. 
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RHO Energy dependent Coulomb phase factor to be multiplied 
into the b coefficients for proton-partial=wave 


mixing. 
A(I,J) The a, b, and c factors occurring in the expression for 
BEL, J) g. LL denotes the coefficient order (order = 21) and 
C(1I,J) J the mixing parameter index. See FIND. 
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Fortran Listing of Program SPECFIT 
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